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O ■ Abstract 

m 

C*~> ' The notion of distance between a global Maxwellian function and 

an arbitrary solution / (with the same total density p at the fixed 
moment t) of Boltzmann equation is introduced. In this way we es- 
sentially generalize the important Kullback-Leibler distance, which 
was used before. Namely, we generalize it for the spatially inhomoge- 
neous case. An extremal problem to find a solution of the Boltzmann 
equation, such that dist{M, /} is minimal in the class of solutions with 
■ the fixed values of energy and of n moments, is solved. The cases of 

the classical and quantum (for Fermi and Bose particles) Boltzmann 
equations are studied and compared. The asymptotics and stability 
of solutions of the Boltzmann equations are also considered. 



1 Introduction. 

We consider the classical and quantum versions of Boltzmann equations 
(where the quantum version contains both the fermion and boson cases). 
The important notion of Kullback-Leibler distance [6], which was fruitfully 
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used before (see further references in the recent papers [UHUED]), ^ s essen- 
tially generalized and new conventional extremal problems, which appear in 
this way, are solved. The solution f(t, x, £) of the Boltzmann equation is 
studied in the bounded domain f2 of the x-space. Such an approach essen- 
tially changes the usual situation, that is, the total energy depends on t and 
the notion of distance between a stationary solution and an arbitrary solution 
of the Boltzmann equation includes the x-space. Thus, the notion of distance 
remains well-defined in the spatially inhomogeneous case too. Recall that the 
Kullback-Leibler distance is defined only in the spatially homogeneous case. 
The comparison of the classical and quantum mechanics, which was treated 
in [T2Tm] . is generalized here for the case of the Boltzmann equations. It is 
especially interesting for the applications that the fermion and boson cases 
are essentially different from this point of view. In the last section of the 
paper we introduce the dissipative and conservative solutions and find the 
conditions under which the stationary solution of the classical Boltzmann 
equation is stable. 

First, we discuss the classical case. The well-known classical Boltzmann 
equation for the monoatomic gas has the form 



where t&R stands for time, x — (xi, . . . , x n )£Q stands for space coordinates, 
C = (Ci, . . . ,£ n )GM n is velocity, and M. denotes the real axis. The collision 
operator Q is defined by the relation 



where B(( — cr)>0 is the collision kernel. Here we used the notation 



= (c* + 0/2 + He* - CI/2, C = (C* + 0/2 - ff|C* - CI/2, (1.3) 



where cxGS' n ~ 1 , that is, aGlR™ and |er| = 1. The solution f(t,x,Q of Boltz- 
mann equation (II. ip is the distribution function of gas. We start with some 
global Maxwellian function M, which is the stationary solution (with the 
total density p) of the Boltzmann equation. The notion of distance between 
the global Maxwellian function and an arbitrary solution / (with the same 
value p of the total density at the fixed moment t) of the Boltzmann equa- 
tion is introduced. As already mentioned before, our approach enables us to 



dt 
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treat also the inhomogeneous case. An extremal problem to find a solution 
of the Boltzmann equation, such that dist{M, /} is minimal in the class of 
solutions with the fixed values of energy and of n moments, is solved. 

The same considerations prove fruitful for the quantum Boltzmann equa- 
tion. Our definition of the quantum entropy S q is slightly different from the 
previous definitions (see [2KTU]). We show that the natural requirement 

Sq^-S c , £—7-0 (S c is the classical entropy) (1.4) 

is not fulfilled in the case of old definition, however ( II. 4p holds in the case of 
our modified definition (see Section E]). 

Some necessary preliminary definitions and results are given in Section |2j 
An important functional, which attains maximum at the global Maxwellian 
function is introduced in Section |3j The distance between solutions of (11. ip 
and the corresponding extremal problem are studied in Section HJ The mod- 
ified Boltzmann equations for Fermi and Bose particles (the quantum cases) 
are considered in Sections |5] and |6j A comparison of the classical and quan- 
tum cases is also conducted in Section [6j Finally, Section [7] is dedicated to 
the asymptotics and stability of solutions. 

We use the notation Cq to denote the class of differentiable functions 
f(C), which tend to zero sufficiently rapidly when ( tends to infinity. 

2 Preliminaries: basic definitions and results 

In this section we present some well-known notions and results connected 
with the Boltzmann equation. The distribution function f(t,x,() is non- 
negative: 

/(£,£, C)>0, (2.1) 

and so the entropy 

S(t, f) = - [ [ fit, x, C) log f(t, x, ()d(dx (2.2) 

is well-defined. 

Definition 2.1 A function <p(() is called a collision invariant if it satisfies 
the relation 

[ 0(C)W,/)(CR = O for all feC 1 ,. (2.3) 
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It is well-known (see that there are the following collision invariants: 
0o(C) = l, UO = d (i = l,2,...,n), n+1 (C) = |C| 2 . (2.4) 

The notions of density p(t, x), total density p(t), mean velocity u(t, x), energy 
E(t,x), and total energy E(t) are introduced via formulas: 

p(t,x)= [ f{t,x,C)d(, p{t)= I p{t,x)dx, (2.5) 



(t,x) = (1/pOM)) / {f(t,x,()d(, (2.6) 



E(t,x)= f ^-f(t,x,()d(, E(t)= f [ ^f(t,x,QdCdx. (2.7) 
The function 

/ = (p/(27rT)"/ 2 ) exp ( - |C - u\ 2 /(2T)). (2.8) 

is called the global Maxwellian and is a function of the mass density p > 
0, bulk velocity u = (ui, . . . ,u n ) and temperature T. We assume that the 
domain Q is bounded and so its volume is bounded too: 

Vol(fi) = V Q < oo. (2.9) 

Therefore, the function 

M(C) = (p/(V a (2nT)^ 2 )) exp ( - |C - u\ 2 /{2T)) (2.10) 
is a global Maxwellian with the constant total density p. 

Proposition 2.1 FTP) / TTie global Maxwellian function M(() is the station- 
ary solution of the Boltzmann equation (11. ip . 

Boltzmann proved in [1] the fundamental result below: 

Theorem 2.1 Le£ f £ Cq be a non-negative solution of equation (II. ID . TTien 
i/ie following inequality holds: 

dS/dt > 0. (2.11) 
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3 Extremal problem 

Similar to the cases considered in [1^115]. an important role is played by the 
functional 

F(f) = {F(f))(t) = XE(t) + S(t), \ = -l/T, (3.1) 

where S and E, respectively, are defined by formulas ( 12. 2 p and (12. 7p . and the 
functional ( 13. ip is considered on the class of functions with the same pit) = p 
at the fixed moment t. The parameters A = — 1/T and p are fixed. 

Now, we use the calculus of variations (see [5|) and find the function / max 
which maximizes the functional ( 13. ip . The corresponding Euler's equation 
takes the form 



S 
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a^-/-/ior/-i /'/ =0. m) 



5f 

Here -jk stands for the functional derivative. Our extremal problem is condi- 
tional and /i is the Lagrange multiplier. Hence, we have 

A^-l-log/ + /i = 0. (3.3) 
From the last relation we obtain 

/ = CV -|C! 2 /(2T)_ (34) 

Formulas ([2~10]) and (IQl imply that 

f = M(Q = , 9 : /0 e~^. (3.5) 

J vs; Vb(27rT)™/ 2 v ; 

In view of (j277|) . and (J3HJ we see that 

F(f)= [ [ L f {t,x,Qd£dz, L f = -(^ + \ogf)f. (3.6) 

For positive / (including the case / = M) and for Lf given in (I3.6p . we have 
the inequality 

Jj&t = - V/ < 0- (3.7) 
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Corollary 3.1 The global Maxwellian function M((), which is defined by 
formula (j3.4jh gives the maximum of the functional F on the class of func- 
tions with the same value p of the total density p{t) at the fixed moment 
t. 

It follows from ([2T5]) . (133]) . and (JEI| that 

F(M » = -" l0g (vW2^)' (3 ' 8 > 

Therefore, Corollary 13.11 can also be proved without using the calculus of 
variation (see [H]). Indeed, taking into account relations (13. 5p . (I3.6p . and 
(13.81) and the fact that the total densities of M and / are equal, we have 

F(M) - F(f) = 1 / M(l - L + llog (3.9) 

Using inequality 1 — x + a; log a; > for x > 0, a; 7^1, we derive from (13.91) 
that 

F(M)-F(/)>0 (f^M). (3.10) 

Remark 3.1 Since the extremal problem is conditional, the connection be- 
tween the energy and entropy can be interpreted in terms of game theory. The 
functional (13. 1 p defines this game. The global Maxwellian function M(() is 
the solution of it. A game interpretation of quantum and classical mechanics 
problems is given in the papers [T4\ , [T5^ . 

Remark 3.2 Inequality (I3.8P is valid for all the non-negative functions f 
with the fixed density p att (not necessarily solutions of the Boltzmann equa- 
tion). 

4 Distance 

Let f(t, x, () be a nonnegative solution of the Boltzmann equation ( 11. ip . We 
assume that T and the value p = p(t) at some moment t are fixed. According 
to (13.101) we have 

F(M)-F(f)>0, (4.1) 
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where the global Maxwellian function M(() is defined in (13. 5p . The equality 
in ( 14.1 p holds if and only if f(t,x,() — Hence, we can introduce the 

following definition of distance between the solution f(t, x, () and the global 
Maxwellian function M(Q: 

dist{M,/} = F(M)-F(/). (4.2) 



Remark 4.1 In the spatially homogeneous case (if not only the total densi- 
ties pM and pf of M and f are equal but the energies Em and Ef are equal 
too ), our definition (14.21) of distance coincides with the Kullback-Leibler dis- 
tance (see \2U$ ). However, our approach enables us to treat also the inhomo- 
geneous case. 

Next, we study the case Em ^ Ef and start with an example. 
Example 4.1 Let T x ^ T and consider the global Maxwellian function 

Direct calculation shows that 

E 1 = E Ml = pnT x /2 ^ E, (4.4) 
F{M X ) = - P flog ( —£— ) - n (l - T l /T)/2] . (4.5) 



v .Vn{2itT x yi* , 

It follows from (jMJ) and (J43]) that 

dist^Mj = —pn[\og(Ti/T) — Ti/T + l)/2. (4.6) 

We introduce the class C(p, E%, U) of non- negative functions functions f(t, x, () 
with the given total density p (see (2.8)), total energy 

f [ ^-f(t,x,C)dCdx = E 1 , (4.7) 

and total moments U — (Ui, U2, ■ ■ ■ , U n ) , where 

Uk = [ [ ( k f(t,xX)d(dx. (4.8) 



Recall that the global Maxwellian function M is defined by ( 13. 5|) . 

Extremal problem. Find a function f , which minimizes the functional 
dist{M, /} on the class C(p, Ei, U). 

The corresponding Euler's equation takes the form 
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(\ + v)\Slf-flogf + l ,f + fJ2~/k(k\ =0. (4.9) 



Recall that our extremal problem is conditional, and p, z/, 7& are the La- 
grange multipliers. Hence, we have 

(A + u) ^ - log / - 1 + fi + 7*C* = 0. (4.10) 

k 

From the last relation we obtain 

/ = C exp ((A + u) ^ + TfcOfc) • ( 4n ) 

k 

According to (12. 5p we have A + v < 0. Now, we rewrite (I4.1ip as 
where 



A + v) V 2 ^ \ \ + u, 

k 



n/2 ST^ 2 

c -= c ( -( A ;, )/2) ./ 2 ^(-^))- < 4 - i3 > 

To calculate the parameters /x, z/, 7*. we use again the well-known formulas 

/DO POO -I 

e~ ae di = ^f^fa, / f e -^ 2 ^= — yf^., a>0. (4.14) 
-oo J— oo ^ a 

Formulas O, (j4T7j), gSJ, f l4TT2|) . and fQ4j) imply that 

C! = p/V n , lk /(\ + v) = -U k /p, -(\ + u)=T l -\ (4.15) 

where 



k 
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Because of (I4.12p and (14.151) we see that / is just another global Maxwellian 
function 

In the same way as (14.51) we obtain: 

F(M2) = ( log (vhp^F 5 ) " " (1 " Tl/r)/2 ) " 2>' 2 ' (4 - 18) 

Moreover, formulas (13. 6p and ( 14. 2 p imply the relations 

dist{M,/}=/ f (L M (t,x,C)-L f (t,x,0)<Kdx, ^-{L M -L f ) = l/f. 

(4.19) 

That is, the functional dist{M, /} attains its minimum on the function / = 
M 2 , which satisfies conditions pit) = p, (14. 7p . and (I4.8p . More precisely, in 
view of (14.181) we have 

1 7T|2 

dist{M,M 2 } = -^(log^/T) - Tx/T + l) + LL. (4.20) 
Hence, the following assertion is valid. 

Proposition 4.1 Let M andM 2 , respectively, be defined by (I3.5P and (14.171) . 

If the function f satisfies conditions pit) = p, (14 ,7p . (14. 8p . and f^M 2 , then 

dist{M, /} > -^(logm/T) - T x /T + 1) + 

Definition 4.1 VFe denote by M the Maxwell function of the form ( 13.51) . 
where p = (l/e)(2irT) n / 2 V n . 

According to (13. 8p we have 

F(M) > F(M), M^M. (4.21) 
Hence the following statement is valid. 
Proposition 4.2 The inequality 

6(f) = F(M) - F(f) > 0, f^M (4.22) 
is fulfilled for all non-negative f . 

We call G in ( 14.22ft the Lyapunov functional, and will study it in greater 
detail in Section [71 
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5 Modified Boltzmann equations for Fermi 
and Bose particles. 

We study the modified Boltzmann equation which takes into account the 
quantum effect 

§£ = -C-V«/ + C(/,/). (5.1) 
The collision operator C is defined by the relation 

C(f, f) = f ! B((- C, a) [f (C) / (CD (1 + e/(C)) (1 + ef (CO) 

- /(C)/ (CO (1 + ef (CO) (1 + ef (C))]dadC (5.2) 

where C' and C* are introduced in ( 11 .31) . and eel. If £ = 0, the right-hand 
side of (15.2p coincides with the right-hand side of (11 .2p . that is, we get the 
classical case. The inequalities e > and e < hold for bosons and fermions, 
respectively. Similar to the classical case the quantum density p e and quan- 
tum energy E £ are given by formulas (12.41) and (12.7ft . respectively. However, 
the quantum entropy S(t,e) (e ^ 0) is defined in a more complicated way: 

S(t,f,e) = - [ [ [f\ogf-(l/e)(l + ef)log(l + ef)+f]d(dx (5.3) 

Remark 5.1 Our definition ( 15. 3ft o/ entropy is slightly different from the 
previous definitions (see f^lTU^). Namely, formula ( 15. 3 j) contains the addi- 
tional summand 

-Pe = ~ I I fd(dx. (5.4) 
We shall show that the natural requirement 

S(e)^S c , e^O (5.5) 
is fulfilled only in the case that (15. 3ft holds. 

6 Modified extremal problem 

1. We assume again that the domain Q is bounded and introduce the func- 
tional 

F £ (f) = XE £ (f) + S(f,e), A = -l/T, (6.1) 
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where E £ (f) and S(f,e) are defined by formulas (12. 7p and (15. 3p respectively. 
The parameters A = —1/T and p are fixed. 

Again we use the calculus of variations (see [2]) and find the function 
fmax which maximizes the functional (16. ip under additional condition 

p(t, x)dx = p. (6.2) 



n 

The corresponding Euler's equation takes the form 

A^- _ i g f + i og (i + £ f) _ i + p = o. (6.3) 

From the last relation we obtain 

f/(l + ef) = Ce-^. (6.4) 

Formula (16. 4p implies that 

ci 2 

Ce 2T 

f = M £ = w . (6.5) 

1 — Cee 2T 

It is required that the distribution M e is positive, that is, 

C > 0, -oo < Ce < 1, (6.6) 

and further we assume that (16. 6p holds. Moreover, (16. 6p yields also the 
positivity of 1 + sM £ : 

M £ (C)>0, l + eM £ (C)>0. (6.7) 

According to (12. 5p and (16.21) . the constant C is defined by the equality 

ci 2 

V a / w dC = p. (6.8) 

jRn 1 _ £ c e -— 

In view of (16. ip . we have the relation which is similar to ( 13 . 6 j) : 

F £ (f)= [ [ L f>e (t,x,()d(dx. (6.9) 
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Though the function L f >£ is more complicated than L f in ( 13. 6 p , we easily get 
an analog of ( 13. 7ft : 

jp L,f = ~~nr+W) K °' (610) 

which clearly holds if / and 1 + ef are positive, including the case that 
f = M £ . 

Corollary 6.1 The functional F £ given by (16. ip attains its maximum (for 
positive functions f satisfying condition (16.81) ) on the function M £ of the form 
(16.51) . That is, for the distance G £ we get 

G £ (f) := F £ (M £ ) - F £ (f) > (f^M £ ). (6.11) 

Remark 6.1 The global Maxwellians M £ play an essential role in boson and 
fermion theories. When the standard approach is used, they appear in a more 
complicated way (see |7[ Ch. V, sections 52, 53] and |3, Ch.l, sections 9, 10]). 

2. Using the spherical coordinates, we calculate the integral on the left-hand 
side of (EHD 

C e ~^ f°° r n - x e'^r 2W 2 

d( = u n -iC I -2-dr, w n _i = (6.12) 



1-eCe-^- ^ l-eCe-2T r W 2 ) 

where a; n _i is the surface area of the {n — l)-sphere of radius 1, and T(z) is 
the Euler's Gamma function. Taking into account (I6.8P and (I6.12p we obtain 

(2nT) n / 2 V n CL n/2 (Ce) = p, (6.13) 

where 



r 2 

oo r n-l e - — 



L " lM = wr^m I — #*' (6 ' 14) 



Because of the equality 



o 



OO 1 

e- ar2 r n ~ x dr = -a- n/2 T(n/2) (6.15) 



the function L n /2(z) admits the expansion 

oo 

L n/2 (z) = Y,{z m - 1 /™ n/2 ), (6.16) 

m=l 

which yields the next statement. 
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Proposition 6.1 The function L n / 2 (z) monotonically increases for 
< z < 1 and 

L 1/2 (l) = L 1 (l) = oo, L n/2 (l)<oo for n > 2. (6.17) 

Remark 6.2 It is easy to see that L n /2(z) = oo for z > 1. 

In view of Proposition 16. II we have: 

Corollary 6.2 If e > (boson case) and either n = 1 or n = 2, then 
equation (16.131) has one and only one solution C such that C > 0, Ce < 1, 
and so (16. 6 j) holds. 

Corollary 6.3 If e > (boson case), n > 2 and 

(27rT) n /V n L„ /2 (l) > ep, (6.18) 

iaen equation (I6.13P aas one and on/?/ one solution C such that C > and 
Ce < 1. J/ ; instead of (16. 18j) . we aawe (27rT) n / 2 Vh£ n /2(l) = ^P; ^en tae so- 
lution of (I6.13P groen fry C = 1/e and the corresponding M e has singularity 
at ( = 0. 

Remark 6.3 The function L n / 2 (z) belongs to the class of the L-functions /2J/ 
and is connected with the famous Riemann zeta-function 

oo 1 
k=l 

by the relation 

L n/2 (1) = C(n/2). (6.20) 
Hence, some useful estimates for L n /2(1) follow. In particular, we get 

£3 /2 (l) = 2.612, L 2 (l) = 1.645, L 5/2 (l) = 1.341, L 3 (l) = 1.202. 

(6.21) 

Let us consider the fermion case (i.e., the case e < 0). The next propo- 
sition easily follows from f !6.14j) and monotonical increase of ax(l + ax)^ 1 
(a > 0) on the positive half-axis. 

Proposition 6.2 Let e < 0. Then the function CL n / 2 (Ce) monotonically 
increases with respect to C > 0. Furthermore, we have CL n / 2 (Ce)— too for 
C-kxi. 

Corollary 6.4 If e < (fermion case), then equation (I6.13P has one and 
only one solution C such that C > 0. 
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3. Consider now the energy for the global Maxwellian M £ : 

f f ICI 2 Ce"w f°° r n+1 e-£ 

E £ (M £ ) = / -™ - s -d(dx/2 = V a v n -iC -dr/2. 

Jn Jm™ i _ E Ce~— Jo 1 - eCe'^r 

(6.22) 

Formulas fl6TT2|) -f l6U4j) and fl6\22|) imply that 

According to ( 14 . 4 j) the corresponding classical energy E — Eq — E c is given 
by the formula 

E C (M) = ^ (M = M ). (6.24) 

Proposition 6.3 If e > (boson case), then we have 

E £ < E c . (6.25) 

If e < (fermion case) and 

either n > 2, — Ce<l or n — 1, -Ce < 3 3/2 /2 5/2 w 0.91, (6.26) 
i/ien we have 

E c < E £ . (6.27) 

Proof. Taking into account (16.161) . we obtain L n / 2+ i(Ce) / L n / 2 {Ce) < 1 for 
e > 0. Hence, in view of (I6.23|) and (I6.24p the inequality (I6.25P holds in the 
boson case. 

If e < and conditions (I6.26P hold, the inequalities 

L n /2{CE) > and -kn/2+l (Ce) - L n/2 (Ce) >0 

follow from f)6.14p and f)6.16p . respectively, and we get 

-W/2+1 (Ce)/L n/2 (Ce) > 1. 

That is, in view of (16.231) and (I6.24p the inequality (I6.27P is proved in the 
fermion case. □ 
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4. For the classical case e = formula (I6.13P (see also ( I2.1(jp ) implies 

C = C = p/V n {2iTT) n/2 . (6.28) 
In view of (pH]l . (ETSj) . and (16T28D we easily derive for M = M that 

S c = ^£ c -plogC . (6.29) 
To calculate the quantum entropy S(M £ , e) we recall (I6.5P and use equalities 
M e = g/(l-eg), 1 + eM £ = (1 - eg)' 1 , g := Ce'^'^ (6.30) 

to simplify the expression, which stands under integral on the right-hand side 
of (15. 3j) and which we denote by Lg: 

L s = M £ (l + log g) + (1/e) log(l - eg). (6.31) 

Substitute log# = logC - (1/2T)|C| 2 into (IQTj) and substitute ( 1PTT) into 
( B to get 

S(M £ ,e) = )-E £ - (1 + \ogC)p - -V n [ log(l - ^R. (6.32) 
Using integration by parts and the definition (12. 7\\ of energy we rewrite f I C3 . 3 2 1) : 
S(M e , e) = ^E £ - (1 + \ogC)p + (6.33) 
From dSHD, (E23), (16T29|) . and fl6T33|) we see that 

77-1-2 

S(M £ , e)-S c = ——{E e - E c ) - plog(CyCo), (6.34) 
ni 

^ - F c = 4^ - s «0 - plog(C/C ) (F c = F ). (6.35) 
nl 

The behavior of C is of interest and we start with the proposition below. 
Proposition 6.4 The following inequalities are valid: 

C>C fore<0; C < C for e > 0. (6.36) 
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Proof. According to (I6.14p and (I6.16P we have 

L n / 2 (zi) < L n/2 (0) = 1 < L n/2 (z 2 ) for z x < < z 2 < 1. (6.37) 

Therefore, it is immediate that 

C L n/2 (C £i) < C , C < C L n/2 (C Q e 2 ) for e x < < e 2 . (6.38) 

In view of Propositions 16. Il and l6.2l the functions CL n / 2 {Cei) and CL n / 2 {Ce 2 ) 
increase with respect to C > 0, and so formulas (16.131) and (16.381) imply 
(IPS]) . □ 

It is immediate from (I6.36P that C is bounded for e > 0. However, C is 
bounded also for the small values of when e is negative. Indeed, let 
-(2C y 1 < e < 0. Then, formula flBTTD yields 

2L n/2 (2C e) > 2L n/2 {-l) > L n/2 (0) = 1. 

Therefore, we have 2CoL n / 2 (2Coe) > Co, which in view of Proposition 16.21 
implies C < 2 Co. 

Now, rewrite (16.131) as z = Cos / L n / 2 {z), where z = Ce, and note that 

4- j c '" £ / s < 1 for \z\ < 1 and small values of e. (Since C is bounded, we see 

;hat |^| < 1/2 for the sufficiently small values of e.) Thus, we apply iteration 
method to the equation z = Cos / L n i 2 {z) and derive 

C = C + O(e), e^O. (6.39) 

Next we note that formula ( 16 . 1 3[) yields CL n / 2 (Ce) = Co- Therefore, taking 
into account (I6.39P we get 

C/Co = l/L n/2 (Ce) = 1 - (Coe)/2 n / 2 + 0(e 2 ). (6.40) 

Moreover, from (I6.40p we see that 

log(C/Co) = -{Coe)/2 n ' 2 + 0(e 2 ). (6.41) 

Using relations f l6U6|l . fl6T23|) . (1Q4) . and fl6T39|l . we derive 

Ee-E c = -^^ + 0(e% e^O. (6.42) 

Because of (16.341) . (I6.35p . (16.411) . and (16.421) . we get the next proposition. 
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Proposition 6.5 For e — > 0, we have equality (I6.42p as well as equalities 
below: 

g(M«, e)-S c = - {n ~ { ^" £ + 0(b\ (6.43) 
I" ^ = 0^ + 0^). (6.44) 

Corollary 6.5 Let S\ < < e 2 be small. Then 

S{M e2 ,e 2 ) < S c < S(M ei ,£i) for n > 2, F £l < F c < F £2 for all n. 

(6.45) 

Remark 6.4 We recall that in view of Proposition fiOl the inequalities 

E £2 <E c <E ei , £l <0<e 2 (6.46) 

hold without the demand for to be small. Here E £2 corresponds to the boson 
and E £l to the fermion case. 

Remark 6.5 We note that relations (16.441) as well as their physical inter- 
pretation are contained in the well-known book by L. Landau and E. Lif- 
shitz Section 55]. 

Conjecture 6.1 Relation (16.271) . which was proved for all —Ce<l (e < 0) 
in the case that n > 2, is valid also for all —Ce<l (e < 0) in the case that 
n = 1. We recall that (I6.27P holds for n = 1 and —Ce < 3 3 / 2 /2 5 / 2 . Moreover, 
(I6.27P holds in the extremal case Ce = — 1. Indeed, using (16.201) . (I6.2ip . the 
relation £(i/2)~ — 1.46 and the well-known equality (see, e.g., |2l P-17]) 

L s (-l) = C( S )(l-2 1 -), (6.47) 

we obtain 

L 3/2 (-l)w0.765, L 1/2 (-l)w0.6. (6.48) 

Hence, L 3 / 2 (—l)/Li/ 2 (—i) > 1 and the conjecture is proved for the case that 
Ce = -1. 
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7 Lyapunov functional 

7.1 Classical case 

In this subsection we prolong the study of the classical Boltzmann equation 
(II. ip and assume that f(t,x,() is its non-negative solution. Using Gauss- 
Ostrogradsky formula we write 

/ / (|C|72)C-V*M^ = / / (\(\ 2 /2)[(-n(x)}fd(da = A(t,n), (7.1) 
Jn Jr™ Jan Jr™ 

I I C-vJd(dx= [ [ [(■n(x)}fd(da = B(t,n), (7.2) 
Jn Jr™ Jan Jr™ 

where dQ is the piecewise smooth boundary of Q, and the integral j dn gda is 
the surface integral with n(x) being the outward unit normal to that surface, 
xedVl. 

Remark 7.1 Here A(t, Q) and B(t, Q) are the total energy flux and the total 
density flux through the surface dQ per unit time, respectively. 

Definition 7.1 We say that a non-negative solution f(t,x,Q of (11.11) be- 
longs to the class of dissipative functions, if A(t, Q)>0 for all t. 

Definition 7.2 We say that a non-negative solution f(t,x,() of (11.11) be- 
longs to the class C(Q) of conservative functions, if A(t, Q) = for all t. 

Clearly we have C(Q)<zD(Q). We note that the same definitions are applicable 
in the quantum case. 

Proposition 7.1 // the inequality f(t,x,() > and condition f(t,x,() = 
f(t,x,—() for xEdQ hold, then we have f(t,x,Q G C(Q). 

Proof. Since / M „(|C| 2 /2)/(t, x, ()(d( = 0, it follows that A(t,Q) = for A 
which is given by (17. ip . □ 

Remark 7.2 The so called bounce-back condition f(t,x,() = f(t,x,—() 
means that particles arriving with a certain velocity to the boundary dQ will 
bounce back with the opposite velocity (see f7P|. p. 16]). 

Corollary 7.1 The global Maxwell functions M of the form (13.51) belong to 
the conservative class C(Q). 
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Furthermore, the next assertion can be easily derived via direct calculation. 

Corollary 7.2 The global Maxwell functions M of the form (I2.10p . also 
belong to C(Q). 

Example 7.1 The well-known and important Maxwellian diffusion example 
(see JX2, p. 16]) is described by the property 

f(t,x,C)=P-{x)M b (0 for xedQ, C-n(x)>0, (7.3) 

where M&(£) has the form (13.51) . When we have 

(\(\ 2 /2)[(-n(x)]fd(dcT> [ [ (\(\ 2 /2)\(-n(x)\fd(da, 

I dil J C_-n{x)>0 J d£l J C,-n(x)<0 

(7.4) 

the function f in (I7.3P is dissipative. If in relation (17.41) we have the equality, 
then f is conservative. Hence, such functions satisfy our statements below 
(and the results below are new even for this case). 

Now, consider the Lyapunov functional G(f) = F(M) — F(f) for the 
equation ( II. ip . According to (I4.22p we have 

G(f) > for f^M, G(M) = 0. (7.5) 
Using Theorem 12.11 we derive the following assertion. 

Theorem 7.1 Let /gCq be a non-negative dissipative solution of (II. ip . 
Then the inequality (dG/dt)<0 is valid. 

Proof. The function <f)(() = \(\ 2 is a collision invariant (i.e., ( 12 .3p holds). 
Therefore, taking into account ( II .ip . ( 17. ip . and Definition 17. II we have 

d 
dt 




(\{\ 2 /2)fd(dx=-A(t,Q)<0, (7.6) 

'ft JR n 



that is, (dE/dt) < 0. Recall also that M is a stationary solution, and so 
dG/dt = —dF(f)/dt. Now, the assertion of the theorem follows from (12.1 ip 
and (I3~TT) . □ 

According to Theorem 17. 11 if its conditions are fulfilled and (G(/))(to) < 
5, then the inequality (G(f))(t) < 5 holds for all t > t . Thus, the following 
important result is proved. 
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Theorem 7.2 If the distance is defined by G and f is dissipative, the sta- 
tionary solution M is locally stable. 

The previous results on local stability [H1EIE] were obtained for the spatially 
homogeneous Boltzmann equation. 



Corollary 7.3 Let conditions of Theorem 7.1 be fulfilled. Then the function 
F(f) monotonically increases with respect to t and is bounded. Hence, there 
is a limit 

lim F(f) = $<F(M). (7.7) 

t— >oo 

Next, assume that the following limits exist: 

Poo = lim pit) ^ 0, U oa = lim U(t), (7.8) 

t— >oo t— >oo 

where pit) and U(t) are given by (12.51) and (j4.8j) . respectively. We see from 
(T3781) and (JESJ that the functions M and M(t) of the form (j3T5) . where p^ 
and p(t), respectively, are substituted in place of p, satisfy relations 

F{M) = log (vsfeO = F(M(t)) - (7 - 9) 



Proposition 7.2 Lei the relations (17. 7ft and ( 17.81) ZioZd. T/ien we /iave i/ie 
inequality 

F(M)-$> |C/ oo | 2 /(2p 0O T). (7.10) 

Moreover, if the inequality (I7.10P torns into equality, there is a unique Maxwell 
function M v of the form (I4.17P (with p = Poo and U = U^) such that 

F(M V ) = $. (7.11) 

If the inequality 1 17. 10ft is strict, that is, 

F(M)-<S>>\U QO \ 2 /(2p QO T), (7.12) 

i/iere are two such functions {Mi and M 2 ) satisfying 

F(M fc ) = $ (A; = 1,2). (7.13) 
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Proof. It is immediate that 

y(x) := x — 1 — logx = for x = 1, y(x) > for x > 0, x ^ 1. 

(7.14) 

Since y > 0, according to Proposition 14.11 we have 

F(M(t)) - F(f(t)) > \U(t)\/(2p(t)T). (7.15) 

In view of (EZD-(E9D and fLUD we get ( 1740]) . 

Now, using (I4.20p we rewrite equation (17. lip (or, correspondingly, (I7.13P ) 
in the form 

F(M)-$-i-^= =x-l-logx, (7.16) 



™Poc V 2 Poo^ 

where or, correspondingly, are expressed via solutions Xk of (I7.16P 
in the form (compare with (I4.17P ): 

T,lf n/ r Poo ( |C - Uoo/P' 

= = iiiwr? exp ( ^jr 

According to (17. 14ft . the equation (17.161) has a unique solution when (IT. 101) 
turns into equality and has two solutions when (I7.10P is a strict inequality. 

□ 



Corollary 7.4 Let the conditions of Proposition pO| be fulfilled. Then 



F{M k ) - F(/)->0, t^oo. (7.17) 



Corollary 7.5 Lei £/ie conditions of Proposition \ 7.£\ be fulfilled, where the 
strict inequality (I7.12p holds. If the limit 

E OD = limE{t) (t-K») (7.18) 

exists and the corresponding solution f(t,x,() converges to a Maxwell func- 
tion, then either = E% or E^ = E 2 . 



Remark 7.3 Proposition 7J? and Corollaries 7^4_ an d ore valid if the 



limit ( 17.7ft exists. We do not suppose there, that the corresponding solution 
f is dissipative. 
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7.2 Quantum case 



Now, we consider the quantum version (15. ip of the Boltzmann equation. The 
corresponding Lyapunov functional G £ (f) has the form (16.111) . 

Theorem 7.3 (see JMj and FTBf ) Let f (/gCq) be a non-negative solution of 
equation (16. ip . Then the inequality 

^>0 (7.19) 

is valid. 

In the same way as in the classical case we obtain the assertions. 

Theorem 7.4 Let / G Cq be a non-negative dissipative solution of equation 
(16.11) . Then the inequality 

<0 (7.20) 



dt 
is valid. 



Corollary 7.6 Let the conditions of Theorem 7.4 be fulfilled. If 

G E (f(t ,x,C))<6, (7.21) 

then 

G e (/(t,z,0)<*, t>t . (7.22) 

Thus, we proved that the stationary solution M £ is locally stable, when the 
distance is defined by G e (f) and the function / is dissipative. 



8 Conclusion 

We see that the study of the Boltzmann equations in a bounded domain Q 
and the suggested new extremal problem allow us to introduce a notion of 
distance and obtain various results for the inhomogeneous classical and quan- 
tum cases. In particular, the notion of the dissipative solutions is introduced 
and asymptotics and stability of solutions of the classical and quantum Boltz- 
mann equations is studied. Following, e.g., [THEE] we pl an al so to consider 
solutions of the Boltzmann equations for the case of Tsallis entropy. The ap- 
proach could be applied to other related equations, such as the Fokker-Planck 
equation. 
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